3  Workshop  “Thermochemical  processes  in  plasma  aerodynamics” 


THEORETICAL  AND  EXPERIMENTAL  INVESTIGATION  OF  A  SHOCK  WAVE 
STRUCTURE  IN  MICROPOROUS  LIQUIDS 


E.A.  Bushmanov,  V.Yu.  Velikodnyi,  V.P.  Vorotilin,  A.V.  Eremeev ,  I.B.  Timofeev*, 

Yu.G.  Yanovsky  and  D.  M.  Van  Wie** 

Institute  of  applied  mechanics  RAS  Moscow,  Russia, 

*M.V.Lomonosov  Moscow  State  University,  Moscow,  Russia,  E-mail:  igortim@orc.ru  , 

** Applied  Physics  Laboratory  The  Johns  Hopkins  University,  USA 

distribution  statement  a 

Approved  for  Public  Release 
Distribution  Unlimited 

Introduction 

Developing  of  new  prospective  samples  of  hypersonic  air  vehicles  with  Mach  number  of 
flight  M>5  requires  a  development  of  new  engines  since  the  efficiency  of  conventional  jet- 
engines  sharply  decreases  in  this  case.  It  is  supposed  to  apply  detonation  and  scram  jets  in  this 
velocity  range.  One  of  the  main  requirement  of  such  engines  developing,  testing  and  running 
consists  in  existence  of  cheap  technological  fuel.  This  fuel  has  to  satisfy  to  ecological 
performance,  it  has  to  ensure  a  stable  work  of  a  scram  jet  at  Mach  numbers  M=3,5-f  14  and  of  a 
detonation  engine  at  M=0^9.  Main  direction  of  investigations  is  connected  with  such  a  fuel 
creation  on  basis  of  the  aviation  kerosene.  In  fact  there  is  a  problem  of  new  fuel  type  creation 
with  controlled  physical  and  chemical  properties  depending  on  flight  conditions. 

It  is  well  known  [1,2]  that  preliminary  bubbling  of  the  liquid  fuel  (for  example  of  the 
aviation  kerosene)  by  an  oxidizer  (air,  oxygen  saturated  air)  or  by  flammable  gases  (CH4 ,  H2, 

C2H2)  essentially  improves  its  spraying  characteristics,  mixing,  range  of  a  stream,  and  it 
decreases  an  ignition  delay  time.  Generating  of  initial  radicals,  products  of  CO,  H2,  H20 
incomplete  burning,  products  of  ZCxHy  conversion  and  cracking  during  this  process  can 
essentially  decrease  an  induction  time  of  chain  reactions  at  lower  initial  temperature  of  the  air- 
fuel  mixture(see.  [2,4]). 

All  these  factors  in  aggregate  favorably  influence  development  of  combustion  and 
detonation  processes.  In  the  work  are  presented  means  and  results  of  new  device  tests,  which 
insure  generating  of  “activated  porous  fuel”  on  liquid  hydrocarbon  basis  and  have  indicated 
above  advantages  as  a  unit. 
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Analysis  of  physical  processes  and  tests  of  devices  for  “activated  porous  fuel” 


generating 

Activated  porous  fuel  generating  process  consists  of  the  following  consequently 
accomplishing  stages: 

1.  Preliminary  bubbling  of  initial  fuel  (for  example  of  the  aviation  kerosene,  straw  oil,  and 
mazut)  under  pressure  of  a  disperser; 

2.  After  transmission  of  it  through  a  cavitator,  where  the  following  bubble  splitting  takes  place; 

3.  Shock  wave  splitting  of  a  gas-disperse  mixture; 

4.  Dispersion  of  this  mixture  to  the  air  stream. 

A  principle  scheme  of  a  device  that  insure  a  realization  of  all  pointed  out  above  processes 
is  represented  in  Fig.l.  All  its  main  units  are  made  of  a  transparent  material-the  organic  glass. 
This  allows  to  undertake  visual  observations  of  processes  inside  this  device  and  with  a  help  of 
these  observations  to  select  optimal  modes  of  its  work.  A  work  of  the  device  takes  place  in  a 
following  way.  An  initial  liquid  (a  fuel)  comes  under  pressure  to  an  area  in  front  of  a  disperser, 
through  which  a  gas  is  forcing  8  (with  a  diameter  of  pores  40h-100  pm)  by  a  tube  10  (a  diameter 
of  a  transmitting  cross  section  is  5h-7  mm).  Created  initial  gas-disperse  stream  comes  to  a 
cavitator  4  (with  the  ratio  of  the  transmitting  cross  section  to  the  width  2,5:1  Oh-4:  15)  through  the 


Fig.l.  The  scheme  of  the  working  part  of  the  set  up:  1  - 
the  aluminum  flange;  2  -  the  tightening  pin;  3  -the 
sprayer  4  -  the  cavitator;  5  -  the  shock  wave 
generatorreHepaTop;  6  -  the  body  of  the  organic  glass; 
7  -  the  metallic  basket;  8  -  the  disperser:  9  -  the  sealing 
device;  10 -the  plastic  tube,  1 1  -  the  transmitting 
working  channel. 
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After  that  the  stream  of  the  gas-disperse  mixture  having  a  velocity  40  -  80  m/s  interacts 
with  a  shock  wave  generator  5  (  a  cylinder  with  the  ratio  of  its  diameter  to  its  length 
3:17-1-5:17),  which  is  located  perpendicular  to  a  flow.  Then  the  gas-disperse  stream  comes  to  a 
sprayer  3  (with  the  ratio  of  its  diameter  to  its  width  2:10-^3.5:10).  A  diameter  of  transmitting 
working  channel  1 1  is  D=18,5  mm. 

The  key  physical  phenomenon,  which  is  used  in  this  device,  is  that  a  sound  velocity  in  a 
gas-disperse  medium  can  be  two  orders  of  magnitude  smaller  than  in  a 
pure  liquid  and  by  an  order  of  magnitude  smaller  than  in  a  gas.  According  to  [5]  the  sound 
velocity  value  in  the  medium  with  the  bubbles,  which  diameter  is  d  ~  0,01  mm,  can  be 
calculated  on  basis  of  the  following  equation: 

2  =  fy=Jc_ _ p _ 

dp  **(1-*)* 

where  a  is  the  sound  velocity,  vis  the  polytrope  exponent,  p  is  pressure,  pis  a  liquid  density, 
(p\  is  a  volume  gas  content.  The  formulae  (1)  is  valid  at  0,1  <  ^  <0,9.  The  sound  velocity  a  at 
<p— >1  works  for  the  sound  velocity  ag  in  the  gas,  but  at  (p  — >  0  it  works  for  the  sound  velocity  in 
the  liquid  auq.  The  minimum  of  the  sound  velocity  is  reached  at  the  volume  gas  content  <p-  0,5. 
For  example,  its  value  is  a  -  23  m/s  at  p  =  1  atm,  and  a=l  m/s  at  p  =  0, 1  atm.  It  appears  from 
this  that  the  stream  velocity  equal  to  the  sound  velocity  is  easily  reached  at  rather  moderate 
pressures  p\o  >5  atm,  and  at  the  volume  gas  contents  at  the  inlet  of  the  cavitator  (p\o  >0,2. 


In  the  photo  2  one  can  see  a  set  up  for  “activated  porous  fuel”  generating  in  working 
condition.  The  shock  wave  front,  appearing  during  the  gas-disperse  stream  flowing  around  the 
cylinder,  is  distinctly  seen.  In  this  case  the  water  was  used  as  a  liquid  and  the  air  as  a  gas.  The 
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gas-disperse  stream  parameters  are  Mi=l,5,  <p,  =  0,7,  Q=40  m/s  (where  Q-  is  a  velocity  of  the 
gas-disperse  stream  CKopocTb). 

It  is  interesting  to  apply  effects  taking  place  at  the  shock  wave  impact  on  the  microporous 
liquid  (fuel)  for  its  preliminary  treatment  by  realization  of  the  conversion  and  cracking 
processes,  generation  of  radicals,  partial  oxidizing  before  spraying  to  a  supersonic  flow  or  a 
detonation  chamber.  Preliminary  estimates  and  calculations  of  processes  in  some  units  of  the 
device  have  been  made  (see  Fig.  1-2).  A  theory  for  a  description  of  a  shock  wave  structure  in  a 
gas-disperse  flow  with  high  volume  gas  content  <p,  =  0,3  +  0,95  has  been  developed.  An 
analytical  solution  of  a  problem  has  been  obtained.  A  theory  of  the  turbulent  microporous  stream 
mixing  with  a  flow  of  the  oxidizer  has  been  developed  as  well. 

A  shock  wave  structure  in  a  microporous  liquid 

The  shock  wave  treatment  of  the  “activated  porous  fuel”  is  the  main  technological  process 
at  its  preparation.  In  the  result  there  can  be  realized  the  isothermal  shocks  (temperatures  of  the 
liquid  and  the  gas  phases  are  equal),  the  adiabatic  shocks  (the  gas  phase  temperature  is 
considerably  higher  than  the  temperature  of  the  liquid  phase)  and  the  intermediate  variants  at  the 
gas-disperse  stream  interaction  with  the  shock  wave  generator  with  respect  to  Mach  number  and 
the  volume  gas  content.  The  model  conceptions  about  the  heterogeneous  medium  as  the  system 
of  the  interpenetrative  components  (phases)  [6]  were  put  to  a  basis  of  the  theory  of  the  shock 
wave  structure  in  the  microporous  liquid.  The  cell  model  of  the  heterogeneous  medium  is  used 
for  the  description  of  the  dissipative  mechanisms.  It  consists  in  the  idea  that  the  definite  volume 
of  the  liquid  “a  cell”  [6]  corresponds  to  each  bubble.  Besides,  the  ideas  and  conceptions 
developed  in  [7,8]  were  used  for  the  description  of  transport  mechanisms  and  dissipative 
processes  in  the  shock  wave.  Their  essence  consists  in  the  following:  separately  each  phase 
represents  two  groups  that  are  characterized  by  their  mean  gasdynamic  parameters.  Each  group 
in  its  turn  consists  of  the  cells  (a  bubble  and  the  volume  of  the  liquid  connected  to  it).  Interaction 
and  the  diffusion  of  the  groups  in  the  shock  wave  front  is  realized  due  to  the  turbulent  transport 
of  the  cells.  The  dissipation  and  energy  exchange  between  the  cells  takes  place  at  the  interface. 
The  additional  condition  used  by  us  at  derivation  of  initial  equations  for  the  description  of  the 
shock  wave  structure  in  the  microporous  medium  is  the  supposition  that  the  enthalpy  of  the 
liquid  does  not  change  in  the  shock  (see.  [5]).  With  account  of  all  said  above  we  have: 


p*(plc,+ps„<pSc„=An, 


(1) 
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Pl(PsC!  Pl(PxC II  ~  ^12'  (2) 

<P!  +<Pr  +<P[+<P[  =  1-  (3) 

pA +pm  +4)+pA,+p,M  +4)=4>  (4) 

A+~Mc'->£f-+^rM 'a =a..  <5> 

2  1  2  /-I 


<*(*/< J  [2(r-D 

4/2  nA^ 


,  [Pitf$ci  +  Pii<$cn  )1 — 


(6) 


1-^ 


v  ci  J 


v^3«fa -^r 

where  An,Al2,Al3,Al4- are  known  constants,  /?;-  is  a  density  of  the  liquid  phase,  pg  ,pgu- are 
densities  of  the  gas  phase  in  the  incident  flow  and  behind  the  shock,  respectively,  c,,c„  -  are  the 
flow  velocities  before  and  behind  the  shock  wave,  P, ,  P„  -  pressure  before  and  behind  the  front, 
<pg,<pg-  volume  contents  of  the  groups  from  the  gas  phas t,<p[,(p[-  volume  contents  of  the 
groups  from  the  liquid  phase-  are  the  variables  dependent  on  the  current  coordinate  -X,  n-  is  the 
polytrope  exponent,  Mr  is  the  Mach  number  in  the  incident  flow,  T,,-  is  know  constant  (see 

[8]),  <p,  -is  the  volume  gas  content  in  the  incident  flow,  ps  =  pg <pg  +  pt(p\\  pT  =  p*n<p\  +  p,(p[- 
So  c,  and/*  can  be  considered  as  constant  equations  (2)-(7)  have  an  analytical  solution.  The 
values  of  c, ,  c„ ,  P, ,  Pu ,  (p, ,  <pn  are  connected  by  the  relation  of  Rankine  and  Hugoniot  type  for 
the  disperse  media  in  this  case  [5,6]: 
j-gr 


<P: 


1  +  exd 


(8) 


lbab  J 


where  X-is  the  current  coordinate,  dbab-is  the  microbubble  diameter, 


lal  [2  nrf(\—(fl) 

I"  11 


i-w" 
<  c>  j 

3 

1  r-3 

1 

T 

:Q 

2(H) 

:,(9) 


-1 

all  the  other  variables  (p[ ,  <pg ,  (pg  can  be  expressed  through  q>[ 

Using  (2)-(9)  one  can  determine  the  pressure  and  the  temperature  distributions  in  the 
microbubbles,  the  density  in  the  shock  wave  front  and  the  shock  wave  width: 
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dbab 

8 


dp 


d(X/dbab)j 


Pn  -Pi 


(10) 


where  8  -  is  the  shock  wave  width.  So  we  have 


dh  „  B 

=  (l-9>/)  ,  (ID 

o  4 

In  Fig.  3  there  is  presented  the  variation  of  the  ratio  of  the  microbubble  diameter  in  the  incident 
flow,  and  the  shock  wave  width  with  respect  to  the  Mach  number  at  different  values  of  the 
volume  gas  content.  One  can  make  conclusions  important  for  practical  applications  on  a  basis  of 
data  presented  in  the  figure.  The  shock  wave  width  is  wide  at  the  gas  content  decreasing 
(p,  <  0,45 ,  and  the  isothermal  shock  is  realized.  The  temperatures  in  the  microbubbles  and  in  the 
liquid  phase  are  approximately  equal  in  this  case.  So  the  fuel  conversion  and  cracking  in  the 
shock  wave  does  not  have  place  even  at  considerable  energy  input  to  a  flow  M ,  —5-5-10. 
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Fig.3.  Dependences  d/6  from  Mach  number  at  different  gas  content,  polytropic  exponent  n 
=  1,4:  1  -  (pi  =  0,9,  2  -  (pi  =  0,8,  3  -  (pi  =  0,7,  4  -  q>i  =  0,6,  5  -  (pi  =  0,5,  6  -  (pi  =  0,4 


The  shock  wave  width  is  small  at  (p,  >  0,7  ,  and  the  adiabatic  mode  of  the  SW  formation  is 
realized  in  this  case.  The  gas  temperature  in  the  bubbles  is  considerably  higher  than  the  liquid 
phase  temperature.  Physical  and  chemical  process  are  intensively  realized  in  the  bubbles  already 
at  the  Mach  numbers  Mi  =4-6.  The  variation  of  the  maximum  achieved  temperature  in  the 
microbubbles  during  the  compression  in  the  shock  wave  with  respect  to  the  Mach  number  is 
presented  in  Table  at  <p,  =  0,8 . 
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Table 


Mi 

2,8 

3,5 

4 

T\\  °K 

962 

1415 

1758 

Estimates  of  the  mixing  parameters  of  the  turbulent  porous  stream  in  the  drifting  air 

flow 

One  the  main  important  technological  processes  during  the  preparation  of  the  “activated 
porous  fuel”  is  its  mixing  with  the  oxidizer  flow.  Preliminary  studies  have  shown  that  the  stream 
spraying  angle  considerably  increases  (up  to  two  times)  at  the  increase  of  the  volume  gas  content 
(p, ,  the  range  of  the  gas-disperse  stream  increases  and  the  thinness  of  the  dispersion  improves 
also.  In  order  to  facilitate  the  search  of  the  model  optimal  working  modes  simultaneously  with 
the  set  up  for  generating  of  “the  activated  porous  fuel”  the  theory  of  the  turbulent  porous  liquid 
stream  interaction  with  the  transversal  gas  flow  has  been  elaborated. 

Formulation  of  a  problem  and  derivation  of  equations  describing  the  turbulent 
stream  motion  in  the  drifting  flow 

The  mechanism  of  the  capture  and  the  derived  from  it  the  capture  velocity  value  for  the 
direct  stream  flow  were  substantiated  earlier  in  [9],  It  was  made  on  a  basis  of  the  conception  of 
the  stream  flow  as  the  continuous  chain  of  the  large  scale  vortices.  At  the  separating  flow  around 
them  the  capture  of  the  liquid  from  the  area  of  the  external  potential  flow  takes  place.  In  the 
present  paper  we  consider  the  generalization  of  the  problem  for  the  case  of  the  turbulent  stream 
flowing  in  the  drifting  flow  of  the  external  gas. 

Let  us  consider  the  following  formulation  of  the  problem.  The  turbulent  stream  of  the 
porous  liquid  flows  out  of  the  circular  hole  of  ro  radius  at  the  initial  velocity  Uo  and  the  density  po 
into  the  space  of  some  external  gas.  This  gas  moves  in  the  direction  of  the  axis  x  with  the  density 
Poo  and  the  velocity  Uoo.  The  principle  scheme  of  the  flow  is  presented  in  Fig.4. 
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Fig.4.  Conventional  scheme  of  the  turbulent 
porous  stream  propagation  in  the  drifted  gas 
flow:ao-the  angel,  at  which  the  stream  comes 
out,  uo  and  po-are  the  initial  values  of  the 
stream  velocity  and  density,  u  and  p-  are  the 
current  values  of  these  parameters,  u»  and  pa>- 
are  the  velosity  and  density  of  the  external  gas. 


The  deviation  of  the  stream  trajectory  takes  place  under  the  action  of  the  external  drifting 
flow  from  its  initial  direction  at  the  angle  oto  to  the  axis  x  with  the  asymptotic  coming  to  the 
direction  along  the  axis  x.  Let  us  introduce  the  curvilinear  coordinate  xj,  which  is  directed  along 
the  tangent  in  each  point  of  the  trajectory 

x.=foV1+y/2 dx> 

where  the  stroke  means  the  derivative  over  the  variable  x. 

Earlier  obtained  in  [9]  equations  of  the  external  media  capture,  of  the  balance  of  the 
longitudinal  momentum  component  balance,  and  the  equation  of  the  stream  composition 
conserve  the  invariant  form  for  the  curvilinear  stream  also.  For  this  purpose  at  the  equations 
derivation  is  necessary  to  use  the  cross  sections  normal  to  the  stream  trajectory  as  it  is  shown  in 
Fig.  4.,  with  the  distinction  that  instead  of  u®  it  is  necessary  to  use  its  projection  uxoo  on  the 
stream  direction  in  the  given  cross  section  xj. 
u  =u  /Jl+y'2 

\oo  oo  V  J 

The  indicated  equations  can  be  written  in  the  following  form,  if  we  take  for  the  typical 
scales  (by  which  the  independent  and  sought  variables  take  the  dimensionless  form)  the  radius  of 
the  initial  stream  cross  section  ro,  its  velocity  uoand  the  density  poin  this  cross  section: 
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d  pur‘ 


dx. 


dx, 


=  2pM]u-uxMjr 
d 


2  2  . 

pu  r  -u 


dx, 


P  —  poo  - 


(Poo-D 


ur 


In  the  obtained  equation  system  we  does  not  have  an  equation  for  the  curvature  of  the 
stream  trajectory.  As  it  was  mentioned  above  the  curvature  of  the  trajectory  takes  place  under  the 
dynamic  head  of  the  external  gas  flow  around  it.  Considering  the  small  part  of  the  stream  as  the 
cylinder  body,  and  making  the  natural  supposition  that  its  streaming  has  the  separating  character, 
we  can  write  the  equation  for  the  resistance  force  acting  a  cylinder  part  of  some  small  length  dL 
in  the  form 

Ff=  yf  p,„u2Xoo(rdL), 

here  yf~  1-  is  the  resistance  coefficient  at  sufficiently  great  Reynolds  numbers.  The  trajectory 
curvature  under  the  action  of  the  resistance  force  leads  to  the  centripetal  force  appearance  that 
counterbalance  it 
Fc  =  pu27t(r2  dL)/R, 

where  R  =  -(1  +  y'2)3/2/y*  is  the  radius  of  curvature  of  the  stream  in  the  considered  point  of  its 
trajectory. 

Equating  them  Ff  and  Fc  to  each  other  one  can  obtain  the  missing  equation  that  closes  the 
system  of  equations,  which  describes  the  stream  motion  in  the  drifting  flow  of  the  external  gas. 


Some  results  of  numerical  calculations 

Obtained  system  of  equations  can  be  solved  in  general  case  only  numerically.  The  solution 
is  dependent  on  three  independent  parameters  -  ratio  of  the  velocity  and  density  of  the  external 
flow  and  the  stream  in  its  initial  cross  section  ux7uo,  pa/po  and  of  the  angle  cio  between  the  initial 
stream  direction  and  the  direction  of  the  external  gas  flow.  The  dependence  of  the  relative 
extension  of  the  stream  part  (at  which  the  stoichiometric  relation  of  the  fuel  composition  and  the 
oxidizer  is  achieved)  of  the  initial  gas  content  (pt  represents  the  most  interest.  Since  the  initial 
gas  content  (p,  in  the  experiments  can  be  varied  in  the  range  from  0,1  to  0,95,  then  for 
Po=(Pk/Pb)(1-^,  )+<p,  we  can  find  the  estimate  p(r=700-r35.  At  the  same  time  the  density  of  the 
mixture  p  differs  from  the  unity  by  the  small  value  about  10'2  at  the  stoichiometric  mass  relation. 
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In  the  Fig.5  one  can  see  the  dependence  of  the  dimensionless  length  Xst  ( q>,  )/X(  (p,  =0)- 
(establishing  of  the  stoichiometric  mixture  composition)  of  the  volume  gas  content  (p,  .  The 
initial  direction  of  the  gas-  disperse  stream  flowing  out  forms  the  right  angle  with  the  direction  of 
the  external  gas  flow. 


0  0.2  0.4  0.6  0.8  s0 

Fig.5  Curves  of  the  dependence  of  stoichiometry 
relation  establishing  length  via  initial  volume  gas 
content  of  the  stream  eo  for  a  number  of  the 
drifting  flow  velocity  values  u<x,,  related  to  its 
value  calculated  at  £o=0  and  u*=l . 

As  it  is  seen  from  the  picture  the  relative  extension  of  the  stream  part  (at  which  the 
stoichiometric  condition  is  achieved)  monotonously  decreases  with  the  increase  of  the  gas 
content  parameter^?, . 

Conclusions 

1.  Means  have  been  elaborated  and  devices  have  been  manufactured  for  production  of 
microporous  liquids  with  the  high  gas  content  and  the  microbubble  diameters  d~  10-50 
|im. 

2.  Shock  waves  visualizing  methods  in  the  microporous  liquids  have  been  elaborated;  photos  of 
the  shock  wave  front  at  the  flowing  around  a  cylinder  have  been  made. 
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3.  The  realizing  possibility  of  the  partial  conversion,  cracking  and  oil  products  oxidizing  in  the 
presented  in  this  paper  devices  on  the  fuel  treatment  has  been  theoretically  validated. 

4.  The  problem  of  the  shock  wave  front  in  the  microporous  liquid  has  been  analytically  solved. 
Modes  of  the  device  work  for  the  generating  of  the  “activated  porous  fuel”  were  forecasted 
with  the  adiabatic  regime  of  the  compression  in  the  shock  were  forecasted  on  its  basis. 

5.  The  theory  of  the  turbulent  porous  liquid  stream  mixing  with  the  transversal  gas  flow 
(oxidizer)  has  been  elaborated.  It  was  shown  on  its  basis  the  possibility  of  the  considerable 
decrease  of  the  combuster  sizes  at  the  increase  of  the  volume  gas  content  in  the  fuel  at  the 
diffusional  mechanism  of  the  combustion. 

References 

1.  Avrashkov  V.,  Baranovsky  S.,  Levin  V.  Gasdynamics  Features  of  Supersonic  Kerosene 
Combustion  //  AIAA-90-5268.  P.8. 

2.  Bushmanov  E.A.,  Velikodnyi  V.Yu.,  Timofeev  I.B.,  Yanovsky  Yu.G.,  David  M.  Van  Wie  On 
the  possibility  of  the  activated  porous  fuel  application  in  processes  of  supersonic  and  detonation 
burning.  International  Symposium  Thermochemical  and  plasma  processes  in  aerodynamics. 
Saint-Petersburg  15-19  July  2002,  Holding  Company  Leninets  P.  8 

3.  Yanovsky  Yu.G.,  Sapgir  G.B.,  Sverdlov  E.D.,  Baikov  A.V.  Comparative  study  of  combustion 
characteristics  of  the  hydrocarbon  fuel  thermal  chemical  decomposition  products  in  the 
submerged  torch.  Symposium  Thermochemical  processes  in  plasma  aerodynamics.  Saint- 
Petersburg  May  30  -  June  3  hjohh  2000.  P.46-52. 

4.  Velikodnyi  V.  Yu.  The  Effect  of  translational  nonequilibrium  in  a  shock  wave  on  initiation  of 
the  H2+CI2  reaction//  Chem.  Phys.  Reports.  2001.  vol.  19(4).  P.  787-792. 

5.  Tsiklauri  G.V.,  Danilin  V.S.,  Seleznev  L.I.  Adiabatic  two-phase  flows.  Moscow.  Atomizdat. 
1973,  444  P. 

6.  Nigmatulin  R.I.  Basis  of  heterogeneous  medium  mechanics.  Moscow.  Nauka.  1978.  336  P. 

7.  Struminskii  V.V.,  Velikodnyi  V.Yu.  A  shock  wave  structure.  Doklady  AS  USSR  .1982. 
V.266.JT9  1.  P.28-31. 

8.  Bashlykov  A.M.,  Velikodnyi  V.Yu.  A  shock  wave  structure  in  a  gas  mixture.  Zhumal 
Technicheskoi  Fiziki.l991.V.  6 1  .Af«  8.  P.  33-42. 

9.  Vorotilin  V.P.  Properties  of  turbulent  streams  obtained  on  the  concept  basis  about  them  as  the 
continuous  chain  of  the  large  scale  vortices.  In  the  collectionof  abstracts  “Stability  of 
homogeneous  and  heterogeneous  liquid  flows”  Ed.  V.V.  Kozlov,  Novosibirsk  2001,  P.44. 


